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Abstract 



> 

OO ■ We prove quantum ergodicity for a family of graphs that are obtained from ergodic one-dimensional 

■ maps of an interval using a procedure introduced by Pakonski a/ (J. P/iys. yl, 34, 9303-9317 (2001)). 

' As observables we take the functions on the interval. The proof is based on the periodic orbit 

expansion of a majorant of the quantum variance. Specifically, given a one-dimensional, Lebesgue- 
measure-preserving map of an interval, we consider an increasingly refined sequence of partitions of 
' the interval. To this sequence we associate a sequence of graphs, whose directed edges correspond to 

■""^ , elements of the partitions and on which the classical dynamics approximates the Perron-Frobenius 

' operator corresponding to the map. We show that, except possibly for subsequences of density 0, the 

eigenstates of the quantum graphs equidistribute in the limit of large graphs. 
fH , For a smaller class of observables we also show that the Egorov property, a correspondence between 

classical and quantum evolution in the semiclassical limit, holds for the quantum graphs in question. 



1 Introduction 

- 1 — I I 

' The quantum ergodicity theorem is one of the central results in quantum chaos. Essentially, it asserts that 

^ I in systems in which the classical dynamics is ergodic the probability measures associated with the squares 

of the moduli of the quantum eigenfunctions converge to the classical invariant measure as one approaches 
the semiclassical limit through almost all sequences of eigenstates (any exceptional subsequences have 
density zero). This was originally proved for flows mElElElEli but it has since been extended to discrete 
dynamical systems (chaotic maps); see, for example [01 [7| |S1 |H1 El (for a very readable introduction to 
the subject, the reader should consult ^J). The methods of proof typically involve applying Egorov-type 
theorems, which relate the time evolution of quantum and classical observables in the semiclassical limit. 

Quantum graphs correspond to associating an operator with a graph. For example, this might be the 
discrete Laplacian acting at the vertices, or the one-dimensional Laplacian acting on functions defined 
on the edges of a (metric) graph, with matching conditions applied at the vertices. Such systems have 
recently been the subject of considerable interest |12) . In particular, quantum graphs have emerged as 
important toy models of quantum chaotic behaviour |13l I14| : if one considers sequences of graphs with 
increasing numbers of edges then, under certain conditions, the quantum eigenvalue statistics converge 
to those of random matrix theory il3_, JJ^ J^S^ JJj JX. _18, 19, _20 . However, relatively little attention has 
been paid to their eigenfunction statistics. For example, quantum ergodicity has not been proved in this 
context. Even though the classical (Markovian) dynamics on a fixed graph is mixing, the difficulty lies 
in dealing with sequences of graphs with increasing numbers of bonds. To date, the only examples that 
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have been studied in this hmit are the star graphs (in which the bonds are connected at a single central 
vertex). However, even though any given star graph is classically ergodic, the limit as the number of 
bonds tends to infinity is not quantum ergodic |^ \2'2i 05| . This is not altogether surprising because the 
star graphs do not satisfy the condition under which one expects the spectral statistics to coincide with 
those of random matrix theory (instead, their spectral statistics coincide with those of integrable systems 
perturbed by a singular scatterer |24[ \2b\ ). It turns out that the star graph eigenfunctions are strongly 
scarred by short periodic orbits (see also |26|). 

The problem of finding examples of sequences of quantum graphs that are quantum ergodic thus 
remains open. It is this problem that we address here. We start by discussing how the question of 
quantum ergodicity on general graphs can be related to the ergodic properties of the eigenvectors of 
an ensemble of unitary matrices. Each ensemble consists of matrices DSq, where So is a fixed unitary 
matrix, determined by the corresponding graph, and D is a random diagonal unitary matrix. 

We then identify a particular sequence of graphs (or matrices Sq) for which quantum ergodicity can 
be established. These are graphs obtained from a construction proposed by Pakonski et al |27j involving 
ergodic one-dimensional maps on an interval. We also prove the analogue of Egorov's theorem for these 
graphs. 

To be explicit, given a one-dimensional, Lebesgue-measure-preserving map S : [0,1] — s- [0,1], we 
consider an increasingly refined sequence of partitions Mn of the interval [0,1]. To this sequence we 
associate a sequence of graphs G„ whose directed edges (bonds) correspond to elements of the partitions. 
The quantum evolution on Gn is described by a unitary matrix U„ such that the corresponding classical 
(Markov) dynamics of G„ approximates the Perron-Frobenius operator associated with S. 

To a classical observable 4> G L^[0, 1] we associate a sequence of quantum observables O„(0) which 
are defined as operators corresponding to multiplication by the average value of (f> on an element of the 
partition. We prove that there is a sequence of sets J„ C {1, . . . , |A^„|} such that 

lim = 1 

and for all sequences {jn}'^=i, jn G Jn, 

Yim^{^pf^ ,0^{<t>)'ipfj) ^ 4>{x)dx, (1) 



where xj}^^^ is the j„-th eigenvector of the graph Gn (compare to the corresponding statement for cat 
maps, |7llllp. This is the analogue of "quantum ergodicity" for the graphs in question. It is equivalent 
to the decay of the quantum variance, 
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in the limit n ^ oo. This equivalence follows from a Chebyshev-type inequality and parallels the textbook 
proof of the statement "for uniformly bounded random variables, convergence in mean square is equivalent 
to convergence in probability" . 

If (j) is Lipschitz continuous, we also prove the Egorov property, 

||U„0„(<^)U;i -O„(0o5)|| ^o(|x„|-i), 

where U„ is the quantum transfer operator corresponding to the graph G„ . The existence of the Egorov 
property provides further justification for the use of the term "quantization" when referring to the se- 
quence U„ obtained from a map S. 

It should be noted that we do not explicitly use the Egorov property (EP) in the proof of quantum 
ergodicity (QE). Even though the more traditional route of deriving QE from EP is available to us, we 
feel that the proof in the present form is likely to be more adaptable to other families of quantum graphs. 
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This paper is organized as follows. In Section 2 we review some of the main issues relating to the 
construction of quantum graphs. In Section 3 we introduce the construction of Pakonski et al |27| and 
proceed to discuss some of its properties. In particular, we prove a sufficient condition for a map to 
be quantizable in the fashion described by |27| . This sufhcient condition, although rather restrictive, 
demonstrates that the class of quantizable maps is sufficiently rich to be interesting. 

In Section 4 we introduce the observables on the quantum graphs obtained from maps. Their quantum 
variance is analyzed in Section 5. In Sections 6 and 7 we prove quantum ergodicity for these observables 
by estimating two different contributions to the variance, and in Section 8 we prove the Egorov property 
in this context. Finally, in Section 9 we discuss some of the issues arising in the proof of these theorems 
and the possibility of extending the proofs to larger classes of graphs. 



2 Quantum graphs 

A quantum graph can be defined in two different, but related, ways. In both constructions we start with a 
graph G — (V, B) where V is a finite set of vertices (or nodes), and B is the set of bonds (or edges). Each 
bond b has a non-zero length, denoted Lh. The lengths Lf, are assumed to be rationally independent. 

The first way to define a quantum graph jl4j is to identify each bond b with the interval [0, Lb] of the 
real line and thus define the L^-space of functions on the graph. Then one can consider the eigenproblem 

-^U6(a;) = X^ub{x). (2) 

This setup has been studied by mathematicians since the 1980s |2H1 1211 1301 1311 and was used in 
physical models prior to that 133 1^ . 

To make the operator in |(2Jl self-adjoint we need to impose matching conditions on the behavior of 
u at the vertices of the graph. One possibility is to impose Kirchhoff conditions:^ we require that u is 
continuous on the vertices, and that the probability current is conserved, i.e. 

V -^Ufc(w) = for all veV, (3) 
^-^ ax 

where the sum is over all bonds that originate from the vertex v (the bonds are now taken to be undirected) 
and the derivatives are taken at the vertex v in the outer direction. The admissible boundary conditions 
were classified in, among other sources, [5^ 157) . 

The second construction considers wave propagation on the graph where each vertex is treated as a 
scatterer and the propagation along the bonds is free. This construction was first considered in ^1] and 
generalized in [^Hl to directed graphs. 

In both constructions one ends up with a unitary matrix S(A) = e*^'^S(0), where L is the diagonal 
matrix of the bond lengths. This matrix gives the eigenvalues {A„} of (0) via the equation 

det(I-S(A„)) -0. (4) 

The dimension of the above matrices is equal to the number B of directed bonds of the graph G. If the 
bonds were initially undirected, each bond is split into two directed bonds of the same length. 

In various sources the notion of the "spectrum a{G) of the graph G" can refer either to the eigen- 
problem (0) (and thus solutions of l@J) or to the eigenphases of the matrix S(A) for an arbitrary A. This 
is not as confusing as it might seem, since the statistical properties of both versions of the spectrum are 
conjectured to coincide when averaged over a large interval of A. 

Similarly, the "eigenvector" of G can refer to one of three objects: 

1. the function u{x) that solves ((2Jl, subject to boundary conditions, for some A G o'(G), 

2. the eigenvector of S(A„) corresponding to the eigenvalue 1, denoted by 0„, 



^sometimes called "Neumann" conditions 
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3. any eigenvector of S(A) for arbitrary A, denoted by '0(A). 

There is a simple correspondence between the first two notions of the eigenvector: the solution u{x) is a 
superposition of plane waves with coefficients given by the elements of Below we discuss a heuristic 
formula which connects the ergodic properties of the second and the third types of eigenvectors. This 
formula provides an additional motivation for the results in the main body of our paper, where we study 
the eigenvectors ■j/'(A). It should be mentioned that these results are fully rigorous and do not rely on 
the heuristic connection. 

To proceed, we need to introduce more notation. By '0fc(A) we will denote the fc-th eigenvector of 
S(A). Our observables are diagonal matrices O acting in the space of directed bonds. The matrix L, as 
before, is the diagonal matrix of the bond lengths. The average bond length, TrL, is denoted by 
L. Quantum ergodicity is the property of almost all eigenvectors to equidistribute. This is equivalent 
to the vanishing of the variance in some limit. For example, we would like to prove that the variance of 
(</>„|0|(/)„) -B^^TrO (and, correspondingly, ('0fe(A)|O|V'fc(A)) -B~^TiO) vanishes. At this point two 
obvious questions arise: (a) with respect to which ensemble is the variance taken, and (b) in which limit 
is it expected to vanish? 

Taking, without loss of generality, Tr O to be zero, we define two variances 

V^iA,B) = =1- ^(0JO|0„)^ 

^ A„<A 

V^iS{X),B) = lf](0,(A)|O|t^,(A))^ 

k=l 

where N{A) — ATrL/27r is the mean number of the eigenvalues in the interval [0, A]. 

A heuristic calculation presented in Appendix El suggests that, if the bond lengths are rationally 
independent. 



iV(A) (0„|L|0„)/L A Jo ^ ' 



(5) 



Thus, if the lengths of the bonds are taken from a narrow distribution, the two variances are intimately 
connected. Moreover, following 39' one can show that the limit on the right-hand side coincides with the 
average of (DS{0), B), where D are uniformly distributed random unitary diagonal matrices. Thus 
equation Q relates the quantum ergodic properties of a graph to the like properties of an ensemble of 
random matrices. 

Equation (|SJ) suggests that one cannot in general expect the variance to vanish in the limit A — s- oo. It 
is natural, however, to expect ergodicity in the limit B ^ oo (cf. 7 ). A serious associated problem here 
is the choice of an appropriate sequence of graphs and observables. One sequence of graphs, the quantum 
star graphs, has been investigated in |21II22| . and it was found, in particular, that the variance V^(A, B) 
does not vanish even when B ^ oo. This is not altogether surprising because the star graphs are known 
to exhibit non-standard spectral statistics |14', '24^ , corresponding to integrable systems perturbed by a 
point-scatterer, rather than to chaotic systems [25| (for a review of the quantum fluctuation statistics of 
star graphs see |^). This is due to the fact that the spectral gap in their Markov transition matrix closes 
more quickly as i? — > oo (like 1/ B) than is the case for graphs exhibiting truly quantum chaotic behaviour. 
The lack of quantum ergodicity for the star graphs is related to the existence of strong scarring of the 
eigenf unctions by short periodic orbits. In the present article we study sequences of graphs generated 
from 1-dimensional maps of an interval in a fashion suggested in |27| . We prove that for a suitable choice 
of observables, the variance {S{X), B) converges to for any A, given that the original 1-diniensional 
map was ergodic. This is a stronger statement than the convergence when averaged with respect to A, 
as suggested by relation ©. 
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3 Quantum graphs obtained from Id maps 



Pakonski et al 27 proposed a procedure to associate a sequence of quantum graphs to a one-dimensional 
map of an interval. In this section we review their construction and proceed to investigate some of its 
properties. 

We consider maps of an interval, which we take to be [0, 1]. A partition M of the interval [0, 1] will 
be taken to mean a finite collection of open disjoint intervals Ej (henceforth called atoms) such that 

i=i 

where M — \M\ denotes the number of intervals in the partition. We will denote by £{A4) the set of 
endpoints of the partition Al. In a slight abuse of the notation we will also denote by A4 the cr-algebra 
generated by the atoms of the partition A4. When considering sequences {A^n} of partitions, each 
partition will be a refinement of the previous one, £{Ain) C £{A4„+i)- We will write A4„ ^ -^n+i to 
describe this statement. 

Condition 1. We consider maps S : [0, 1] [0, 1] that satisfy the following conditions: 

(a) the Lebesgue (uniform) measure fi is preserved by the map S: /i(A) = ii (S'^^(j4)) for any measur- 
able set A; 

(b) there exists a partition Mq of the interval [0, 1] into Mq equal atoms, with S linear on each atom; 

(c) the set of endpoints £{Ai) is forward-invariant under the action of S: S{£) C £. 

An example of a map S satisfying Conditionals shown on Figure ^ The tent map with slope 2 is 
another such example. 

Remark 1. The conditions on the map S imply, in particular, that for any two atoms E and E\ either 
S{E) is disjoint with E' or S{E) D E' . 

Remark 2. It is possible to generalize the construction to maps that preserve a measure different from 
Lebesgue, but such a map would have to be topologically conjugate to a piecewise linear map satisfying 
the above properties. To maintain a degree of generality we strive to make explicit the conditions that 
are imposed on the map 5* and the measure /i. 

The Frobenius-Perron operator, reduced to measures constant on each atom of a partition M., can be 
described by a matrix B of size M = The entries of the matrix are given by 

j,{E,nS-HE,)) 

- — Wi) — ' 

and can be described as the answer to the question "what proportion of the set Ej gets mapped into 
-Efc" . An example of a map S and the corresponding matrices B for two different partitions are shown 
on Figure n 

If we view the interval [0, 1] with the uniform measure as a probability space, we can write Bjk = 

¥{S{x) € Ek\x G Ej). 

Lemma 1. Let the set of endpoints of the partition M be invariant under S. Then the matrix B defined 
by ^ satisfies the following properties: 

1. B is stochastic, 

M 
k=l 
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Figure 1: An example of a quantizable map and the corresponding matrices B: (a) \A4\ ~ 4, the atoms 
of the partition are represented by the vertices of the graph; (b) \A4 \ ~ 8, the atoms of the partition are 
represented by the edges of the graph. 



2. If the atoms Ej of the partition M. have equal measure and the map S preserves this measure then 
B is doubly stochastic 

M 

3. If the atoms Ej of the partition M. have equal measure and if the map S is linear with respect to fi 
on each atom Ej (i.e. fi{S{A)) — C^{A) for some C and all A € Ej) then 



_l^{ft-^,S-^{E,^)) 

■'jlj2 ■ ' ' ^jk-ljk ~ 

Proof. Part n follows directly from © 



M(£^j-n[o,i]) 



1. 



Partialis similar: if fJ.{Ej) — m for all j then 

M 1 / / A/ 



i=i \ \j=i 



Part 121 is a consequence of the fact that, if ^{S{A)) — Cj^{A) for all A C Ej, then Bjk is ei- 
ther or l/Cj. Consider first the case Bj^j^^^ = for some r. By definition of B, this means that 
fi {Ej^ n S-\Ej^^,)) = 0. Therefore, 

f,{s-^{E,^)ns-^-'{E,,.^,))=o, 
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and the expression on the right hand side of Q evaluates to zero. 

Now consider the case n — 2 (the case of general n being analogous) with both Bj^j-^ and Bj^j^ being 
different from zero. Let, to simplify the notation, jr — r. Then 

^I{Eons-^ {E.ns-^jE^))) 
fi{Eins-HE2)) 

where the last equality is true by virtue of linearity of S. Using the definition of Bi 2 and the identity 
S-^ (AnB) = S-^ (A) n S-^ (B) , we arrive to 

f,{Eons-\E,)ns~\E2)) 

which is the sought result, given that n{Ei) = /i(£^o)- D 

As mentioned earlier, we are interested in sequences of partitions. 
Condition 2. We consider sequences of partitions A^„ that satisfy 

(a) the atoms within each of the partition have equal measure; 

(b) the sets of endpoints £{Mn) are forward-invariant under the action of S; 

(c) the set of the endpoints of A4„ contains the j-th pre-image of the endpoints of Aio for all j = 1, . . . n. 

Remark 3. Given a map S satisfying Condition ^ one can always construct a sequence of partitions 
satisfying Condition [21 

Remark 4. Conditions [ ]|^b)| and mc)| imply that the map S is non-contracting, ^{S{A)) > fJ.{A). If the 
map is ergodic (see definition ^ in sectional). Condition Ua)| implies that the prcimages of £{A4o) with 
respect to S are dense in [0, 1]. This, in turn, implies that the size of the atoms of the partitions A^„ 
tends to zero (or, equivalently, Af„ 00). 

The following lemma explains the way in which such sequences of partitions 'resolve' the dynamics. 

Lemma 2. Given a partition Al„ satisfying Condition\^ let kg and fc„ be such that 5"(£^fc„) D E^^^ (cf. 
Remark^. Then there exists a unique sequence fci, . . . , kn-i such that 

X € Ek„ and S"{x) e Ek„ => S^x) € Ek^ VO < j < n. 



Proof. Consider an atom E of the partition Mn- Condition [^c]] means that for every j = 0, . . . fc, the 
image (E) lies in a single atom of the "primary" partition Ma. Since the map S is one-to-one on each 
atom of Ada, we conclude, by induction, that is one-to-one on E for every j = 1, . . . n -f 1. 

Assume that the statement of the lemma is incorrect: there are two points, x and y, that satisfy, 
without loss of generality, x,y € Ei, S^{x) e E2, S^{y) e E3 and S^{x), S"-{y) S E4. Remark implies 
the following inclusions: 

S'''{Ei)DE2, S\E,)dE3, S^-^{E2)dE^, S"-'^{E3)d E^. 

Thus each z E4 has n — r-preimages in both sets E2 and E3 and, therefore, two distinct n-preimages 
in El . This contradicts the earlier conclusion that S" is one-to-one on Ei . □ 

Remark 5. Obviously, Lemma |21 is valid if, instead of the "position" of S'^{x) (i.e. the atom Ej such 
that S'"(x) € Ej), we know the position of S'^{x) for some m < n: we can still recover positions of all 
iterates S^{x) for < j < m. In fact, a careful inspection of the proof reveals that the Lemma would 
still be true for m = n + 1. However, if we know only that x G Ei and S'^~^'^{x) G E2, we would not be 
able to pinpoint S^{x), < j < n + 2, to any particular atom of the partition Mn- 
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The next lemma exhibits the block structure of the matrix B. 

Lemma 3. For a partition Mn, n > Q, define an equivalence relation between atoms by setting Ej ~ Ej. 
if S(Ej) intersects S{Ek) and then completing by transitivity. Then the maximum number of elements in 
an equivalence class is uniformly bounded with respect to n. 

For example, in the partition of Figure^] part (b), the atoms Ei, E2, £"3 and E5 form one equivalence 
class and the other four atoms form another equivalence class. Note that, if the atoms of a partition are 
represented by edges of the graph, the equivalence classes correspond to the groups of edges ending in 
the same vertex. For the map in Figure Q the uniform bound on the size of an equivalence class is 4, as 
will be evident from the proof. 

Proof. Take an atom A of the primary partition Mq and let {xi, yi), . . . , (xk, yk) be the disjoint intervals 
forming the pre-image of A with respect to S. By Condition Q Hc) | these intervals contain no endpoints of 
Alo, therefore the map 5* is linear on each interval. Condition [T] also implies that all slopes of the map 
5" are integer. Denote the slope of S on the interval {xj,yj) by Sj. To simplify the notation we assume 
that all Sj are positive. Let p be the least common multiple of sj. 

Condition Ql^c) | implies that Xj and yj are endpoints of the partition A^„. Choose x'^ € £{A4n) such 
that the interval (xj,x'j) contains exactly p/sj atoms of the partition Mn- Since the atoms of A4„ have 
equal length (which we denote by fin), 

Sji {x'j, - Xj, ) = Sj.^ {x'j^ ~ Xj^ ) = PM« , for any ji , ■ (8) 

Moreover, the selected points x'j are the closest to the respective Xj to satisfy both condition ^ and 
x'j e £{Mn)- In particular, this imphes that x'j < yj, since setting x'j = yj would also satisfy condition 

From the above we can conclude that S maps all intervals {xj,x'j) to the same subinterval of A. The 
atoms of Mn making up the intervals {xj,x'^) thus form an equivalence class of size p/si + ■ ■ ■ + p/s^, 
which is independent of n. 

We can now repeat this procedure with intervals {x[,yi), . . . , {x'f.,yk) and, thereafter, with all atoms 
A of the partition Mq. If some of the slopes Sj are negative, the procedure would still go through with 
minor variations. □ 

Having obtained a sequence of doubly stochastic matrices B„ we define their "quantizations" as 
unitary matrices U„ such that 

iB„),k = \iUnhkf. (9) 

The doubly stochastic matrices B for which finding a corresponding U is possible are called unistochastic. 

Condition 3. We assume that the map S is such that all of the corresponding matrices B„, bar finitely 
many, are unistochastic. 

Not all bistochastic matrices are unistochastic. However, formulating a general sufficient conditions 
that ensure unistochasticity is a question of considerable difficulty. The interested reader is referred to 
[271 14(J| and the references therein where some necessary conditions are discussed and where examples 
of maps satisfying and failing Condition |31 are given. To convince the reader that the class of maps S 
satisfying Condition |21 is far from empty we state the following sufficient condition. 

Lemma 4. If the slopes of the map S satisfying Condition^ are all equal (modulo sign), Condition\^is 
also satisfied. 

Proof. This Lemma follows simply from the proof of Lemma |31 Indeed, let s be the absolute value of the 
slope of S. Then all matrices B„, n > 0, have a block structure with blocks of the size s x s and elements 
1/s. Thus the question is really about finding an s x s unitary matrix with all elements satisfying l/Tjfcp = 
1/s. One example of such matrix is the Fourier matrix with elements Ujk — exp{27rijfc/s}/Y^. □ 

Example 1. An example of a map S which has unequal slopes but is still unistochastic is provided by 
the map of Figure ^ 
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Remark 6. An observant reader would notice that, given one unitary U satisfying @, one can produce 
infinitely many such matrices. For example, one can multiply a given U by an arbitrary diagonal unitary 
matrix. However, the results of our paper do not depend on the precise choice of matrices U„, provided 
that condition Q is satisfied. 

One can associate a graph to the matrices B and U in the following way: the indices of the matrices 
enumerate the directed edges of the graph; the end of an edge j coincides with the start of the edge k 
if the matrix element Bjk is non-zero. The number of distinct vertices in such a construction should be 
maximized, then the vertices will correspond to the equivalence classes of Lemma 01 

The matrix B defines a Markov chain on the edges of the graph with Bjk representing the transition 
probability from j to k. The matrix U can be viewed as a quantum propagator on the graph. This 
geometrical interpretation of the two matrices as a graph will be helpful in the later sections when we 
use trajectories on the graph to describe properties of the eigenvectors of U. 

It is also possible to associate vertices of a graph to the indices of B, see Figure ^ part (a). We use 
directed edges for reasons of tradition, rather than convenience. 



4 Quantization of the observables 

Having defined the sequences of unitary matrices U„, ergodic properties of whose eigenvectors we are 
going to study, we need a final ingredient, the observables 0„. For a general sequence of graphs, it is 
not obvious how to define a consistent sequence of observables. In our case, however, there is a natural 
answer. 

We use the discretizations of functions (j) E L^[0, 1] as our observables. Fix a partition M (the semi- 
classical limit corresponds to \A4\ oo). If the function </> is constant on each atom of the partition 
A^, its quantization O = O(0) is a diagonal matrix with entries Ojj = (^{x) where x £ Ej. If ^ is not 
constant on the atoms of A^, we replace (j) by its local average. More precisely, we introduce the piecewise 
constant function (j) defined by 



?(^) = /i N / Hy)dfJ'{y), where Ej 3 x. 
K^j) Je, 



Then we define O — O(0) as before, by 

On - O«(0) = 4>m = ^ md^^{y). (10) 

It is convenient to describe (j) using the notions of probability theory. In probabilistic language, </> is a 
random variable defined on the probability space [0, 1] and (p is its conditional expectation, = E [0| M]. 
We will also use the notation of expectation to denote the integral over [0, 1]: 

Ecl)^ [ <j>{x)dn{x). 



In particular, \\(f>\\2 = {M^f^^- 

To prove quantum ergodicity, we will rely on the ergodicity of the classical map S. Since our ob- 
servables are in L^, the relevant version of the ergodic theorem is the L"^ ergodic theorem (see, e.g., 

mi)- 

Definition 1. A map S : [0, 1] [0, 1] is ergodic if any set A C [0, 1] satisfying 

S-\A) = A 

has either full or zero measure. 
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Theorem 1. (L^ Ergodic Theorem) If 4> E L^[0,1\ and S is ergodic, then 

](l}o ~eA ^0. (11) 



\ t=i 



Since the ergodic theorem appUes to any function from S L^, it also apphes to (j), whatever partition 
A4 was used to produce it. Unfortunately, a uniform estimate for the rate of convergence in (|ll|l for 
different hat- versions of the same (p is not known 03 ■ However, it is easy to see that, for fixed T, 

Vt(0) ^ VtW (12) 

as the partition in the definition of </> gets finer. 



5 Quantum variance of an observable 

Given a map S and a sequence of partitions Mn we have constructed a sequence of Markov matrices B„, 
which, in turn, give rise to unitary matrices U„. On the other hand we are given an observable and 
we have constructed a corresponding sequence of diagonal matrices 0„, which "quantize" (j). We denote 
by Mn the number of atoms in the partition Mn- This is also the size of the matrices B„, U„ and 0„. 
The semiclassical limit corresponds to M„ — > oo. 

Let where j = 1,...,M„ denote the orthonormal eigenvectors of U„. If an eigenvalue is 

degenerate, the particular choice of the basis for its eigenspace is unimportant. As discussed earlier, to 
show quantum ergodicity it is sufficient to prove that 



and 



V- = ^t.\ H"'' - E(<^)|' ^ 0, (14) 



3 = i 

as n — > OO. It is straightforward to verify (|13|l . Indeed, from the unitarity of U„ and the definition of 



1 / si 

i / , (n) „ , (n)\ i 



Mn ^ 



0„t/>5")) = _ Tr 0„ = E (^) = E(0). 



Thus the main task is to show (|14|l . 

Without loss of generality we can assume that E(0) = 0. In what follows wc will omit the sub- and 
super-scripts n unless we want to underline the dependence of a quantity on n and on the partition Al„. 

To obtain an estimate of Vn we employ some standard manipulations. If xp is an eigenvector of a 
unitary matrix U, we have, for any matrix O (not necessarily diagonal) and all t G N 

(ip, Oip) = (uV, ouV) = (u*)*ou*v) • 

Summing this equality over t = 0, . . . T — 1 we obtain 

We introduce the shorthand On,T for the time average of 0„, 

1 

0".T = ^5](u:)*o„u:,. 

t=o 
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Using Cauchy-Schwarz inequality and ortlionormality of {ipj} wc estimate 

|(V,OV')|' = |(V,OtV')|' < (OtV,OtV) = (V,OJOtV), 



and obtain 



= ]i7-5:i(^.'OV',) 



^ (V„ O^Ott/',) = ^ TV (OJOt) '=^' Kin, T) 



(15) 



It is important to note that the above inequahty is vahd for all values of T. Thus, to show that 
— > 0, we are free to choose an appropriate T = T{n) for each n as long as we can demonstrate that 



0, 



as n — *■ oo. In the following sections we prove that T{n) = n is a suitable choice for this task. 
For our purposes, it is more convenient to work with the matrices S„,t defined by 



T-l 



t=0 



which is equivalent to working with Or since SJ^St = (U^Ot)*(U^Ot) = O^Ot- 
Multiplying S^-S out we obtain 



1 1 ^ 

-Tr(S^ST) = ^ E \{St)sJ?. 



M 



We can expand the entries of St in terms of trajectories on the graph. Using the definition of St, we 
obtain 



T-l 



(St)./ 



t=0 bo,-bT 

^ ^ Ubo,bi ■ ■ ■ UbT._i,bT I ^ ^ ^ Obj.-t,bT-t 1 ) 



6n,...5T 



*=0 



where the inner sum in the first line is over all sequences of bonds satisfying bo = s and ^ f ■ Such a 
sequence of bonds we will call a trajectory. Only trajectories compatible with the graph's geometry (i.e. 
those for which Ubjbj+i 7^ 0) contribute to K{n,T). A trajectory r = (60, ... , br) is said to have length 
T and amplitude 

= Ubobi ' ' ' ^bT-ibT • 

We will denote by the average of the observable over the trajectory r, 



I dci 1 

*r = ^ (Obi 61 



ObT^T ) ' 



To summarize, we have shown that 



1 1 
F„ < if (n, = Tr (SJSt) = E 



,/=i 



r:s— >/ 



]t^E E (16) 
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where the inner sum is over all possible trajectories of length T starting at s and finishing at /. 



6 Diagonal terms 

Equation H16|l is reminiscent of a trace formula expansion of the spectral form factor (i.e. of the Fourier 
transform of the spectral two-point correlation function)^, in particular of a graph, see e.g. ,14. . Such 
expansions are notoriously difficult to analyze rigorously as both T and the size of the graph increase. 
The starting point of any such analysis is the evaluation of the contribution from the diagonal terms ^ 
obtained by restricting the last sum in (|16|l to identical trajectories, ti = T2. It is usually assumed that 
the off-diagonal terms sum up to a subdominant contribution, when T and the size of the graph scale 
appropriately. This idea, called the diagonal approximation was first introduced for a general class of 
systems in 03]. On graphs it was explored, in particular, in ^^^]. It is difficult, however, to give an 
a priori estimate on the size of the off-diagonal contributions and the analysis is usually restricted to 
evaluating the contributions coming from specific classes of interacting trajectories jl5[ll6[[T7lll8j . 

Our strategy now is to calculate the contribution from the diagonal terms in (|16|l . Then we will show 
that, in the case of graphs constructed from Id maps, we can actually estimate the off-diagonal terms by 
virtue of being able to choose an appropriate T = T(n). 

To evaluate the diagonal contribution 



we make two observations. First, by the definition of the amplitude A^- and the defining property of the 
matrix U, equation lO, we obtain 

I^T-I^ = \Ubo,bi P • • • \UbT-i.,bT P -^bo.fcl ' ' ' Bbj.-l,bT- 

Now we recall Lemma ^ partial and conclude that 

, _i^{rk=oS-\Eb,)) 

On the other hand, by definition of 

T T T 

<f r = ^ E 0.„.. = ^ E ^{Eb, ) - ^ E O 5* (5-* (i?, J) , 
t=l t=l t=l 

where <j>{Eb) denotes the (constant) value of the function cf) on the atom Eb. In fact, it is easy to see that 
i'x coincides with the value of the function 

on the set H^Lo '^~*(-^''t) £'f)o,- -,6T i if this set is non-empty. If it is empty, the value of <I>^ is of no 
consequence since the trajectory r is then incompatible with the graph's geometry and At = 0. 

The measure of all atoms Eb is assumed to be equal. More precisely, it is equal to 1/Af„, since iV/„ is 
the total number of the atoms. Collecting our observations together, we can express the diagonal term 
as 

,2 \l(Eba,...,hT) 



Vt(0). 



^To underline this similarity we have used in 1151 the traditional notation for the form factor, K 
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Thus, by the ergodic theorem (Theorem^, K'^'^^^^\n,T) goes to zero as T — > oo. On the other hand, 
K{n,T) is bounded below by a non-negative Vn which is, generically, non-zero for a fixed n. This shows 
that the diagonal term is a poor approximation to K{n,T) in the limit T — > oo. Luckily, this is not the 
limit we have to take. 



7 Completion of the proof of quantum ergodicity 

Lemma 121 has a very important consequence for the inner sum in IjlBf) . 

Lemma 5. The diagonal term X^'^'^e) (n, T) gives the exact value of K{n,T) up to time T — n, i.e. 

X(n,T) = i^('^''^«)(n,r) = Vt(?) if T < n, 

where = E [(/)| Mn]- 

Proof. By Lemma |21 for every pair of bonds s and /, there is at most one trajectory going from s to f 
in T < n steps. Thus, for T <n, 



K{n,T) = — y 



sj=l 



J2 '^rAr 
r:s^f 



is I*' 



As a consequence, we have the following result. 



□ 



Theorem 2. (Quantum Ergodicity) Let the map S and the sequence of partitions {Mn} satisfy 
Conditions^ anrf[3 let {U„} be the corresponding sequence of unitary matrices with eigenvectors 
■j/'^"''; and let {0„} be a sequence of diagonal matrices corresponding to an observable (f> G L^[0, 1] with 
E0 = 0. If S is ergodic, then 



Vn = 



as n oo. 



Proof. The variance Vn is majorized by K{n,T) for any T. We combine Lemma|Slwith equation (|12|l we 
conclude that, for a fixed T, 

K{n, T) Vrict)) as n ^ oo. 

Now we use the standard e/2 argument: for any e > 0, by Theoremnwe can find T such that Vt(</>) < e/2. 
Having fixed this T, we find n{e,T) such that \K{n,T) - Vt(0)| < e/2 for aU n > n{e,T). Combining 
the above, 

Vn < K{n,T) < e/2 + e/2 

as long as n > n{e,T). Since e was arbitrary, we conclude that Vn ^ 0. □ 

Remark 7. One can avoid the e/2 argument in the following way. Taking the limit n oo of the 
inequality Vn < K[n,T) produces 

< limsupK < limsup i^(n,T) = Vt(<?!'). 

n — >00 71 — >OQ 

Now taking the T oc limit, we obtain 



< limsupy^ ~ lim sup lim sup < limsupVT(0) — 0. 



13 



8 Egorov property 

In Section^lwe defined a procedure to obtain a piecewise constant function given a function (f> G L'^[Q, 1] 
It is enlightening to see how (f)o S is related to 0. 
By definition, 



(j) o S 



Since 5 is linear on Ej, its derivative is constant. In fact, it is easy to see that — Bjk, where 

y € Ej and S{y) G Ek- Thus we have 



(p o S 



where the sum is over the decomposition of the set S{Ej) into atoms Ek- Since Bjk = whenever 
Ek n S{Ej) is empty and since n{Ej) is independent of j, we arrive to the following conclusion 

Lemma 6. //, for a given partition M., the matrices B, 0{(j)) and 0((/>o 5) are defined according to 0) 
and 1^1U\) then 

M 

0.j.j{cj)oS) = y^-Bj-fcOfcfc('^), 

k=l 

where M is the number of atoms in the partition M. . 

Lemma is a rather beautiful manifestation of the inter-consistency between the discretization pro- 
cedures for maps S and observables (p G L'^. Namely, the discretization commutes with the action of S 
on L^. In this. Lemma |H| is a classical analogue of the Egorov property, a result which shows that the 
unitary matrices U„ faithfully represent the action of the classical map S. 

Theorem 3. (Egorov property) Let the map S and the sequence of partitions {Mn} satisfy Conditions 
^\^and\^ let {U„} be the corresponding sequence of unitary matrices with eigenvectors ip'j^'' ; and let 
{On} be a sequence of diagonal matrices corresponding to an observable (p. If 4> is Lipschitz continuous 
then 

||U„O„(0)U-i - O„(0o 5)11 = 0(Af„-'), 
where the norm is the operator norm on the Euclidean space R*^" . 

Proof. We fix the partition A^„, denote the corresponding UOU^^ by Q and observe that, while 0{(f)oS) 
is a diagonal matrix, Q is not necessarily so. 

First we treat the diagonal elements of Q. Writing them out explicitly we get 

Af„ Af„ A/„ 

r—1 r—1 r—1 

where we used the unitarity of U and its defining property, |t/j>P = Bjj. and Lemma El 
For the off-diagonal elements of Q we have 

A/„ M„ A/„ 



r=0 r=0 
A/„ 



r=0 r=0 
A/„ 



r=0 
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where C is any constant and we have used the unitarity of U to conclude that the second sum is zero. 
We estimate, using Cauchy-Schwarz, 



M„ 

IQjfel < max \Orr - C*! ^ \UjrUkr\ 

< maxlOrr ^ C\ — max \(j)(x) — C\. 

r xeS{E,)nS{Ek) 



If (j) is Lipschitz continuous with 



dcf \(t>{x) - (j){y)\ 
\up = sup < oo, 



iLip 



we can estimate further, by choosing appropriate C, 

\Qjk\ < ^\m^,^fi{S{E,) n S{Ek)) < i||0||Lipmax|5'(x)|A^(i?,) (x . 

Since Ujr is non-zero only if S{Ej) n 7^ and Ukr is non-zero only if S{Ek) D Er ^ 9, we conclude 
that Qjk = if S{Ej) and S{Ek) are disjoint. Thus the matrix Q is of block-diagonal structure, each 
block corresponding to an equivalence class as defined by Lemma |31 The norm of Q — O(0 o S") is equal 
to the maximum of the norms of the blocks. A norm of a block, in turn, is bounded by its dimension 
times the maximum absolute value of the element of the block. The dimension of a block is uniformly 
bounded by Lemma |3| Thus we get 

||Q-O(0o5)||<i?(5)M[idP 

for some constant -0(5*) which is independent of and n. □ 
Remark 8. If the function is only assumed to be continuous on [0, 1], one can prove a weaker property: 

||U„0„(,/))U-i - O„(0o ^)|| ^0 as M„ ^ 00. 



9 Discussion 

We have succeeded in proving quantum ergodicity (QE) for a special class of sequences of quantum 
graphs. However, we would like to mention that the result is expected to hold for much broader class of 
graphs. 

It is true that, given a finite quantum graph G, one can associate a Id map to it by reversing the 
process described in the paper. Thereafter, it is possible to produce a sequence of graphs, one of which 
will coincide with the original graph G, and answer the question of QE for this sequence. In this sense, 
each graph corresponds to a Id map. However, this is not true for every sequence of graphs. In fact, it is 
not true for most sequences. Examples of such sequences include star graphs with Kirchhoff conditions at 
the central vertex (for which the question of QE has been answered negatively), the complete (Kirchhoff) 
graphs, and the star graphs with Fourier central vertex for both of which the QE is expected (but 
is not known) to hold in some form. 

It is reassuring that the proof of QE in the present article suggests a direction for possible gener- 
alizations: study the diagonal terms and then find an estimate for the off-diagonal ones. However, for 
the sequences of graphs described above the diagonal approximation ceases to be exact for T > 1 (cf. 
Lemma O. This makes estimation of the off-diagonal terms a much more difficult task. 

Another interesting question to consider is whether quantum unique ergodicity (when the convergence 
in happens along all sequences of eigenvectors) is true for any quantum graphs. This has been answered 
in the negative 26 for graphs with Kirchhoff vertices but is unclear for other types if graphs. 
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A Connection between variances and 

To demonstrate relation (jSJ we start with summarizing the notation introduced in Section [3 Let the 
unitary B x B matrix S be defined by S = S(A) = e*''''^S(0), where L is the diagonal matrix of the 
bond lengths of the graph and S(0) is some fixed unitary matrix. Let {A„} be the (real) solutions of the 
equation det(I — S(A)) = 0. We assume that the spectrum {An} is non-degenerate, which is a generic 
situation |^. 

Denote by ^„ the normahzed eigenvector of S(A„) corresponding to the eigenvalue 1. By ■0fc(A) we 
denote the k-th normalized eigenvector of S(A). We further denote by e*^'=^^-' the eigenvalues of S(A), 
with 9k chosen to be continuous (indeed smooth) functions of A. 

When A — A„ there is an index k for which ^^(A) = mod 2tt. For this index k we also have 

= V'fe(A). 

Let A be a self-adjoint matrix (a generalization of the observable O) with trace (without loss of 
generality). We are interested in the relationship between two variances. 



and 

B 



y^(S(A),i3) = 



B 

fc=i 

where N{A) — ATrL/27r is the mean number of the eigenvalues A„ in the interval [0, A] 
Introducing the notation 

we observe that 



Tr AS™(A) = ^^fee™''''-(^l (17) 



fc=i 



In particular, Tr A = X^^^i ^fc — 0. From the theory of distributions we know that 



\ m=l 



converges, in the limit £ — > 0, to 

oo 

m= E S{e-27Tr), 

r— — oo 

where 5 is the Dirac delta function. Substituting in the above identity 9 = 9k, multiplying by Ak and 
performing the summation over k yields 



J2^eiek)Ak ^^(f2 e^"" (Tr AS'" + Tr AS-'") J . 

fe=l \m=l / 



(18) 
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As e ^ this converges to 

where, given A„, k is chosen to satisfy Ok ~ 0. It is shown in Jl] that 

Clearly, 0^(A„) > and so we can drop the modulus around ^fc(A„) in the previous equation. Now we 
need the following properties of the approximants to the Dirac delta function, 

27r lim eAiix) = A(x) 



and, if xi X2, 
Applying these identities gives 



27r lim e {A^{x — xi) + A^{x — X2) f — A(x — xi) + A{x — X2)- 



2 



27r Ihn e J2 AM^k - ^ = E TTT)^^^ " ^"^^ 

\k=l I k=l n=l ''^ 

Integrating the right-hand side with respect to A we get 

We will use this quantity to approximate (A, B) . It is a good approximation if the bond lengths of 
the graph are approximately 1 (i.e. the matrix L is approximately unity): 

L^iS^iKB) < V^{A,B) < L,„axV^^(A,B), 

where £max and imin sue the maximal and minimal bond lengths. 
Using p8|l and expanding the square we obtain 

V^(A,B) = =J— [ 2TT\ime AJek)Ak] dX 
N{A) Jo -0 

= —= / lime V e-^(l™il+l"^l)TrAS™i TrAS^^dA. 



mi ,m2 — — ctt 



We now take the A — > oo limit of the above expression and interchange it with the e-limit. Due to the 
rational independence of the bond lengths, the limit 

1 

lim = / TrAS™iTrAS"-^dA 



A^oo N{A) Jo 

is zero whenever mi ^ — m2. Indeed, we recall that S — S(A) — e''^'^S(0) and expand the trace 
Tr AS™ = J2 ^b„MS{0)b,M ■ ■ ■ S{0)b„^M exp(zALp), 

6o,...,b™ 

where Lp — L^-^ + • • • + Lb^^ . If Lp is a sum of mi terms and Lq is a sum of TO2 "t-i terms, they cannot 
be equal. Thus, the only case when the phase factors exp(iALp) in a product of two traces can cancel 
each other is when mi = —m2- 
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We arrive at 

rA 



lim V^(A. B) = Km ^ / Km e V e-^"*^ Tr AS™ Tr AS"'" d\ 

L^oo ' A^oo 7riV(A) Jo ^1 

1 ^ 
= lim ^ / y Al{\)d\ 

A— 27riV(A)7o 



fe=i 

f-A 



A^cc> 27riV(A) 
Here we used expansion (|17ll and the fact that 



Km / y^(S(A),i?)dA. 

7riV(A) Jo 



Kme y e-2™^e^"'(*'=i-^''2) = hm 

m— 1 

is 1/2 when Ok^ = Ok^ and otherwise. 
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